Composition-tuned smeared phase transitions 
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Phase transitions in random systems are smeared if individual spatial regions can order indepen- 
dently of the bulk system. In this paper, we study such smeared phase transitions (both classical 
and quantum) in substitutional alloys Ai-^B^ that can be tuned from an ordered phase at compo- 
sition x = to a disordered phase at x — 1. We show that the ordered phase develops a pronounced 
tail that extends over all compositions x < 1. Using optimal fluctuation theory, we derive the com- 
position dependence of the order parameter and other quantities in the tail of the smeared phase 
transition. We also compare our results to computer simulations of a toy model, and we discuss 
experiments. 
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I. INTRODUCTION 

When a phase transition occurs in a randomly disor- 
dered system, one of the most basic questions to ask is 
whether the transition is still sharp, i.e., associated with 
a singularity in the free energy. Naively, one might expect 
that random disorder rounds or smears any critical point 
because different spatial regions undergo the transition 
at different values of the control parameter. This expec- 
tation turns out to be mistaken, as classical (thermal) 
continuous phase transitions generically remain sharp in 
the presence of weak randomness. The reason is that a 
finite-size region cannot undergo a true phase transition 
at any nonzero temperature because its partition func- 
tion must be analytic. Thus, true static long-range order 
can only be established via a collective phenomenon in 
the entire system (see, e.g., Ref. Q] for a pedagogical dis- 
cussion). 

Recent work has established, however, that some phase 
transitions are indeed smeared by random disorder. This 
can happen at zero-temperature quantum phase transi- 
tions when the order parameter fluctuations are over- 
damped because they are coupled to an (infinite) heat 
bathi^ As the damping hampers the dynamics, suffi- 
ciently large but finite-size regions can undergo the phase 
transition independently from the bulk system. Once 
several such regions have developed static order, their lo- 
cal order parameters can be aligned by an infinitesimally 
small mutual interaction. Thus, global order develops 
gradually, and the global phase transition is smeared. 
Classical thermal phase transitions can also be smeared 
provided the disorder is perfectly correlated in at least 
two dimensions. In these cases, individual "slabs" of 
finite thickness undergo the phase transition indepen- 
dently of the bulk system^ 

The existing theoretical work on smeared phase tran- 
sitions focuses on situations in which a sample with 
some fixed degree of randomness is tuned through the 
transition by changing the temperature (for classical 
transitions) or the appropriate quantum control param- 
eter such as pressure or magnetic field (for quantum 
phase transitions). However, many experiments are per- 



formed on substitutional alloys such as CePdi-^Rhj, or 
Sri-a-Ca^RuOs. These materials can be tuned from an 
ordered phase (ferromagnetic for the two examples) at 
composition x = to a disordered phase at x — 1 while 
keeping the temperature and other external parameters 
fixed, i.e., they undergo a phase transition as a function 
of composition. The composition parameter x actually 
plays a dual role in these transitions. On the one hand, 
x is the control parameter of the phase transition. On 
the other hand, changing x also changes the degree of 
randomness. If such a composition-tuned phase transi- 
tion is smeared, its behavior can therefore be expected 
to be different than that of smeared transitions occurring 
at fixed randomness. 

In this paper, we investigate the properties of 
composition-tuned smeared phase transitions in substi- 
tutional alloys of the type Ai-^B^. We show that the 
ordered phase extends over the entire composition range 
x < 1, and we derive the behavior of the system in the 
tail of the smeared transition. Our paper is organized 
as follows. In Sec. [Til wc consider a smeared quantum 
phase transition in an itinerant magnet. We use optimal 
fluctuation theory to derive the composition dependence 
of the order parameter, the phase boundary, and other 
quantities. In Section ITU we briefly discuss how the the- 
ory is modified for smeared classical transitions in sys- 
tems with correlated disorder. Section HVl is devoted to 
computer simulations of a toy model that illustrate and 
confirm our theory. We conclude in Sec. [V] by comparing 
composition-tuned smeared transitions with those occur- 
ring at fixed randomness. We also discuss experiments. 



II. SMEARED QUANTUM PHASE 
TRANSITION 

A. Model and phase diagram 

In this section we investigate the ferromagnetic or anti- 
ferromagnetic quantum phase transition of itinerant elec- 
trons with Ising order parameter symmetry. In the ab- 
sence of quenched randomness, the Landau-Ginzburg- 
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Wilson free energy functional of this transition in d space 
dimensions reads^i 

S = J dydz il>(y)T(y, z)i}j{z) + u J dy ^(y) . (1) 

Here, ip is a scalar order parameter field, y = (y, r) com- 
prises imaginary time r and <i-dimensional spatial posi- 
tion y, J dy = J dy Jq dr, and u is the standard quar- 
tic coefficient. T(y, z) denotes the bare inverse propaga- 
tor (two-point vertex) whose Fourier transform reads 

r(q,w„) = r + Coq 2 +7o(q)|w„| . (2) 

Here, r is the distance from criticality^ £o is a micro- 
scopic length scale, and uj n is a Matsubara frequency. 
The dynamical part of T(q, w„) is proportional to |cj„|. 
This overdamped dynamics reflects the Ohmic dissipa- 
tion caused by the coupling between the order param- 
eter fluctuations and the gapless fermionic excitations 
in an itinerant system. The damping coefficient 7o(q) 
is q-independent for an antiferromagnetic transition but 
proportional to l/|q| or l/|q| 2 for ballistic and diffusive 
ferromagnets, respectively. 

We now consider two materials A and B. Substance A 
is in the magnetic phase, implying a negative distance 
from criticality, r A < 0, while substance B is nonmag- 
netic with rs > 0. By randomly substituting B-atoms 
for the A-atoms to form a binary alloy Ai^B^, we can 
drive the system through a composition-driven magnetic 
quantum phase transition. 

A crucial role in this transition is played by rare A- 
rich spatial regions. They can be locally in the magnetic 
phase even if the bulk system is nonmagnetic. In the 
presence of Ohmic dissipation, the low-energy physics of 
each such region is equivalent to that of a dissipative two- 
level system which is known to undergo, with increasing 
dissipation strength, a phase transition from a fluctuating 
to a localized phased Therefore, the quantum dynamics 
of sufficiently large rare regions completely freezes^ and 
they behave as classical superspins. At zero temperature, 
these classical superspins can be aligned by an infinites- 
imally small residual interaction which is always present 
as they are coupled via the fluctuations of the paramag- 
netic bulk system. The order parameter is thus spatially 
very inhomogeneous, but its average is nonzero for any 
x < 1 implying that the global quantum phase transi- 
tion is smeared by the disorder inherent in the random 
positions of the A and B atomsi 2 ' 11 ! 12 

At small but nonzero temperatures, the static mag- 
netic order on the rare regions is destroyed, and a finite 
interaction of the order of the temperature is necessary 
to align them. This restores a sharp phase transition 
at some transition temperature T c (x) which rapidly de- 
creases with increasing x but reaches zero only at x = 1. 
If the temperature is raised above T c , the locally ordered 
rare regions act as independent classical moments, lead- 
ing to super-paramagnetic behavior. A sketch of the re- 
sulting phase diagram is shown in Fig. [TJ 
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exponentially power-law tail 

FIG. 1. (Color online) Schematic temperature-composition 
phase diagram of a binary alloy K\- x Yi x displaying a smeared 
quantum phase transition. In the tail of the magnetic phase, 
which stretches all the way to x = 1, the rare regions are 
aligned. Above T c , they act as independent classical moments, 
resulting in super-paramagnetic (PM) behavior. x° c marks 
the critical composition in average potential approximation 
defined in (0. 

B. Optimal fluctuation theory 

In this section, we use optimal fluctuation theor y 13 i 14 
to derive the properties of the tail of the smeared quan- 
tum phase transition. This is the composition range 
where a few rare regions have developed static magnetic 
order but their density is so small that they are very 
weakly coupled. 

A crude estimate of the transition point in the binary 
alloy Ai_ x B x can be obtained by simply averaging the 
distance from criticality, r av = (I — x)r A -\-xrB- The tran- 
sition point corresponds to r av = 0. This gives the critical 
composition in "average potential approximation," 

x° c = -r A /(r B - r A ) . (3) 

Let us now consider a single A-rich rare region of lin- 
ear size Lrr embedded in a nonmagnetic bulk sample. 
If the concentration x\ oc of B atoms in this region is be- 
low some critical concentration x c {Lrr), the region will 
develop local magnetic order. The value of the critical 
concentration follows straightforwardly from finite-size 
scaling,^^ 

x c {Lrr) = x° c - DL-% , (4) 

where 4> is the finite-size shift exponent and D is a con- 
stant. Within mean-field theory (which should be qual- 
itatively correct in our case because the clean transition 
is above its upper critical dimension 6 ), one finds tfi = 2 
and D — £ 2 /(rs — r A )- Since x c (Lrr) must be positive, 
((4]) implies that a rare region needs to be larger than 
£min = [D jx^) 1 !^ to develop local magnetic order. 

As the last ingredient of our optimal fluctuation theory, 
we now analyze the random distribution of the atoms in 
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the sample. For simplicity, we assume that the lattice 
sites are occupied independently by either A or B atoms 
with probabilities 1 — x and x, respectively. Modifications 
due to deviations from a pure random distribution (i.e., 
clustering) will be discussed in the concluding section fVl 
The probability of finding Nb = Nx\ oc sites occupied by 
B atoms in a spatial region with a total of N 
is given by the binomial distribution 



l rr sites 



P{N,x loc ) = 



N 
N B 



(1-x) 



N-Nb x N b 



(5) 



We are interested in the regime x > x° c where the bulk 
system will not be magnetically ordered but xi oc = 
Nb/N < x c (Lrr) such that local order is possible in 
the region considered. 

To estimate the total zero-temperature order param- 
eter M in the tail of the smeared transition (where the 
rare regions are very weakly coupled) , we can simply sum 
over all rare regions displaying local order 



M 



dL 



RR 



dx ioc m(N, xioc) P(N, xioc) 



(6) 

Here, m(N, x\ oc ) is the order parameter of a single region 
of N sites and local composition xi oc ; and we have sup- 
pressed a combinatorial pref actor. We now analyze this 
integral in two parameter regions, (i) the regime where 
x is somewhat larger than x° c but not by too much, and 
(ii) the far tail of the transition at x — Y 1. 

If x is not much larger than x® , the rare regions are ex- 
pected to be large, and we can approximate the binomial 
distribution ([3]) by a Gaussian, 



P(N, Xioc) = 



y/2nN{l - N) 



exp 



-N 



2x(l 



(7) 



To exponential accuracy in x, the integral (j6|) can now be 
easily performed in saddle point approximation. Neglect- 
ing m(N, xioc), which only modifies power-law prefactors, 
we find that large rare regions of size L RR = D(2<fi — 
d)/[d(x — x®)] and maximum possible B-concentration 
xioc = x c ~ DL RR dominate the integral. Inserting these 
saddle point values into the integrand yields the compo- 
sition dependence of the order parameter asil 



M ~ exp 



-C 



(x - x») 2 ~ rf /^ 
x(l — x) 



(8) 



is a non-universal 



where C = 2(D/d) d /*(20 - df^ 
constant. 

Let us now analyze the far tail of the smeared tran- 
sition, x — > 1. In this regime, the binomial distribution 
cannot be approximated by a Gaussian. Nonetheless, the 
integral (|6]) can be estimated in saddle-point approxima- 
tion. We find that for x — >■ 1, the integral is dominated 
by pure-A regions of the minimum size that permits local 



magnetic order. This means Lrr = L Tl 



(D/x 



OU/0 



and xioc = 0. Inserting these values into the integrand 
of we find that the leading composition dependence 
of the order parameter in the limit x — > 1 is given by a 
non-universal power law, 



M ~ (1 - x) 1 



(1 



){D/x°) d/ <t' 



(9) 



We thus find that M is nonzero in the entire composition 
range < x < 1, illustrating the notion of a smeared 
quantum phase transition. 

So far, we have focused on the zero-temperature order 
parameter. Other quantities can be found in an analo- 
gous manner. Let us, for example, determine the phase 
boundary, i.e., the composition dependence of the critical 
temperature T c . As was discussed in Sec. Ill Al the static 
magnetism of the rare regions is destroyed at nonzero 
temperatures. Consequently, magnetic long-range order 
in the sample can only develop, if the rare regions are 
coupled by an interaction of the order of the tempera- 
ture. The typical distance between neighboring locally 
ordered rare regions can be estimated from their den- 
sity, p, as r typ ~ ~ Af _1 / d . Within the Landau- 
Ginzburg- Wilson theory (|H2p . the interaction between 
two rare regions drops off exponentially with their dis- 
tance r, -Ei n t ~ exp(— r/£b), where £{, is the bulk corre- 
lation length. This leads to a double-exponential depen- 
dence of T c on x for compositions somewhat above x®, 
i.e., ln(l/T c ) ~ exp{C(x - x°) 2 ~ d ^/[dx(l - x)]}. For 
x — > 1, we find ln(l/T c ) ~ (1 — x)~ L ™n/ d . However, in 
a real metallic magnet, the locally ordered rare regions 
are coupled by an RKKY-type interaction that decays 
as a power law with distance, E ln t ~ r" d , rather than 
exponentially — (This interaction is not contained in the 
long- wavelength expansion implied in ©.) Therefore, 
the composition dependence of the critical temperature 
takes the same form as that of the magnetization, 



T c ~ exp 



-C 



(x - x° c ) 2 - d /<t> 



x(l — x) 

for compositions somewhat above x° and 

\££.„ _ ft ~\(D/x a c ) d/,i > 



T, 



(1 - x)^» = (1 - x) ( 



(10) 



(11) 



in the far tail of the smeared transition, x — > 1. 

We now turn to the order parameter susceptibility. It 
consists of two different contributions, one from the para- 
magnetic bulk system and one from the locally ordered 
rare regions. The bulk system provides a finite, non- 
critical background throughout the tail of the smeared 
transition. Let us discuss the rare region contribution 
in more detail. At zero temperature, the total order pa- 
rameter M is nonzero for all x < 1. The rare regions 
therefore always feel a symmetry-breaking effective field 
which cuts off any possible divergence of their susceptibil- 
ities. We conclude that the zero-temperature susceptibil- 
ity does not diverge anywhere in the tail of the smeared 
transition. If the temperature is raised above T c , the rel- 
ative alignment of the rare regions is lost, and they be- 
have as independent large (classical) moments, leading 
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to a super-paramagnetic temperature dependence of the 
susceptibility, x ~ 1/T (see Fig.Q}. At even higher tem- 
peratures, when the damping of the quantum dynamics 
becomes unimportant, we expect the usual non- universal 
quantum Griffiths power- laws, \ ~ T 1 ) where A is the 
Griffiths exponent . 11 ' 12 ! 19 



for x — > 1. The same substitution of d by d± was also 
found for smeared classical transitions tuned by temper- 
ature rather than composition.— 

IV. COMPUTER SIMULATIONS 



III. SMEARED CLASSICAL PHASE 
TRANSITION 

Classical (thermal) phase transitions with uncorrelated 
disorder cannot be smeared because all rare regions are of 
finite size and can thus not undergo a true phase transi- 
tion at any nonzero temperature. However, perfect disor- 
der correlations in one or more dimensions lead to rare re- 
gions that are infinitely extended in the thermodynamic 
limit. If the number of correlated dimensions is high 
enough, these infinitely large rare regions can undergo 
the phase transition independently of the bulk system, 
leading to a smearing of the global phase transition. 4 
This happens, for example, in a randomly layered Ising 
magnet, i.e., an Ising model with disorder correlated in 
two dimensions. 5 

In this section, we discuss how the theory of Sec. 
HTl is modified for these smeared classical phase transi- 
tions. For dcfinitcness, we consider a classical Landau- 
Ginzburg- Wilson free energy in d dimensions, 

S = J dy ^(y)[r - d*]rjj(y) + u J dy ^(y) . (12) 

As in the quantum case, we now consider a binary "al- 
loy" Ai_ x B x of two materials A and B. The atoms are 
arranged randomly in d± dimensions, while they are per- 
fectly correlated in d|| = d — dj_ dimensions. For example, 
if d± = 1 and d|| = 2, the system would consist of a ran- 
dom sequence of layers, each made up of only A atoms 
or only B atoms. 

If the correlated dimension d\\ is sufficiently large, the 
"alloy" undergoes a smeared classical phase transition as 
the composition x is tuned from to 1 at a (fixed) tem- 
perature at which material A is magnetically ordered, 
ta < 0, while material B is in the nonmagnetic phase, 
t\b > 0. The optimal fluctuation theory for the behavior 
in the tail of the smeared transition can be developed 
along the same lines as the theory in Sec. [TTJ The only 
important difference stems from the fact that the ran- 
domness is restricted to d± dimensions. The dimension- 
ality d in eqs. ((5]) and © therefore needs to be replaced 
by d±, leading to 



M ~ exp 



-C 



(X - £0)2-^/0 



x(l — x) 

for compositions somewhat above x° c and 

M ~ (1 - x) L & = (1 - x)^/^* 



(13) 



(14) 



To verify the predictions of the optimal fluctuation 
theory in Sec. |TT] and to illustrate our results, we have 
performed computer simulations of a toy model, viz., a 
classical Ising model with d space-like dimensions and 
one time-like dimension. The interactions are between 
nearest neighbors in the space-like directions but infinite- 
ranged in the time-like ones. This (d+l)-dimensional toy 
model retains the possibility of static order on the rare 
regions (which is crucial for the transition being smeared) 
but permits system sizes large enough to study exponen- 
tially rare events. The Hamiltonian reads 



H 



— Sy yT S Z . T ' ~ JySy. T Sy^ T ' (15) 



(y,z>, t,t' 



y,i~,T' 



Here y and z are d-dimensional space-like coordinates 
and t is the time-like coordinate. L T is the system 
size in time direction and (y, z) denotes pairs of near- 
est neighbors on the hyper-cubic lattice in space. J y is a 
quenched random variable having the binary distribution 
P(J) = (1 - x ) S{J - J h ) + x 5{J - Ji) with J h > Ji. 
In this classical model L T plays the role of the inverse 
temperature in the corresponding quantum system and 
the classical temperature plays the role of the quantum 
tuning parameter. Because the interaction is infinite- 
ranged in time, the time-like dimension can be treated 
in mean-field theory. For L T — > oo, this leads to a set of 
coupled mean-field equations for the local magnetizations 



m y = 0-/L T ) J2t Sy.r- They read 

niy — tanh/3 [J y m y + i 



h] 



(16) 



where the sum is over all nearest neighbors of site y and 
h — > is a very small symmetry-breaking magnetic field 
which we typically set to 10~ 12 . If all J y = Jh, the system 
undergoes a (sharp) phase transition at Th = Jh + 2d, 
and if all J y = J;, it undergoes the transition at Ti — 
Ji + 2d. In the temperature range Th > T > TJ, the 
phase transition can therefore be tuned by composition 
x. 

The mean-field equations (HHJ) can be solved efficiently 
in a self-consistency cycle. Using this approach, we stud- 
ied systems in one, two, and three space dimensions. The 
system sizes were up to L— 10000 in Id, and up to L = 100 
in 2d and 3d. For each parameter set, the data were aver- 
aged over a large number of disorder realizations. Details 
will be given with the individual results below. 

Figure [2] shows an overview over the magnetization M 
as a function of composition x for a (3 + l)-dimensional 
system at several values of the classical temperature in 
the interval Th > T > T/. The figure clearly demon- 
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FIG. 2. (Color online) Magnetization M vs composition x for 
a (3 + 1) -dimensional system having Jh = 20, Ji = 8 and sev- 
eral values of the classical temperature T. The data represent 
averages over 100 samples of size L — 100. The values of the 
critical concentration in "average potential approximation," 
shown for comparison. 



strates that the magnetic phase extends significantly be- 
yond the "average potential" value x° c = (T/j — T)/{Th — 
Ti). In this sense, the magnetic phase in our binary al- 
loy benefits from the randomness. In agreement with the 
smeared phase transition scenario, the data also show 
that M(x) develops a pronounced tail towards x = 1. 
(By comparing different system sizes, we can exclude that 
the tail is due to simple finite-size rounding^) We per- 
formed similar simulations for systems in one and two 
space dimensions, with analogous results. 

To verify the theoretical predictions of the optimal fluc- 
tuation theory developed in Sec. [Til we now analyze the 
tail of the smeared phase transition in more detail. Fig- 
ure [3] shows a semi-logarithmic plot of the magnetization 
M vs. the composition x for a (l + l)-dimcnsional system, 
a (2 + l)-dimensional system, and a (3 + l)-dimensional 
one. In all examples, the data follow the theoretical pre- 
diction ((SJ over at least 2 orders of magnitude in M in a 
transient regime of intermediate compositions x. 

We also check the behavior of the magnetization for 
compositions very close to x = 1. Since © predicts a 
non-universal power law, we plot log(M) vs. log(l— x) for 
a (3 + l)-dimensional system in Fig. @] The figure shows 
that the magnetization tail indeed decays as a power of 
(1 — x) with x — > 1. The exponent increases with in- 
creasing temperature in agreement with the prediction 
that it measures the minimum size iV m i n ~ -C'min a rare 
regions needs to have to undergo the transition indepen- 
dently. The inset of Fig. 2] shows a fit of the exponent to 
L d min ~ K{T)]-^ = \{Th—T)/ (T/,,— T;)]~ 3 / 2 . Theequa- 
tion describes the data reasonably well; the deviations at 
small exponents can be explained by the fact that our 
theory assumes the rare-region size to be a continuous 
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FIG. 3. (Color online) log(M) vs x in the tail of the transition 
for three example systems: (i) (3+l)-dimensional system with 
L = 100, J h = 20, Ji = 8, and T = 23, (ii) (2 + l)-dimensional 
system with L — 100, Jh = 15, Ji = 8, and T = 18, and (hi) 
(1 + l)-dimensional system with L = 10000, Jh = 11, Ji = 
8, and T = 12.8. All data are averages over 100 disorder 
configurations. The solid lines are fits to ((8|, with the fit 
intervals restricted to x £ (0.25,0.55) in (1+1) dimensions, 
(0.6,0.72) in (2+1) dimensions and (0.7,0.82) for the (3+1)- 
dimensional example. 

variable which is not fulfilled for rare regions consisting 
of just a few lattice sites. 

Our computer simulation thus confirm the theoretical 
predictions in both composition regions in the tail of the 
transition. In a transient regime above x®, we observe 
the exponential dependence ||5J) while the magnetization 
for x — > 1 follows the non-universal power law (|9|). 



V. CONCLUSIONS 

In summary, we have investigated phase transitions 
that are tuned by changing the composition x in a ran- 
dom binary alloy Ai^^B^ where pure A is in the ordered 
phase while pure B is in the disordered phase. If in- 
dividual, rare A-rich spatial regions develop true static 
order, they can be aligned by an infinitesimal residual 
interaction. This results in the smearing of the global 
phase transition^in agreement with the classification put 
forward in Ref. 120 . 

As an example, we have studied the quantum phase 
transition of an itinerant Ising magnet of the type 
Ai-ajB^. At zero temperature, the ordered phase in this 
binary alloy extends over the entire composition range 
x < 1, illustrating the notion of a smeared quantum 
phase transition. Upon raising the temperature, a sharp 
phase transition is restored, but the transition tempera- 
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iog 10 (i-x) 

FIG. 4. (Color online) log(M) vs log(l - x) for a (3 + 1)- 
dimensional system with L — 100, Jh = 20, Ji — 8 and several 
temperatures. All data are averages over 100 disorder config- 
urations. The solid lines are fits to the power-law ([9]). The 
inset shows the exponent as a function of temperature, with 
the solid line being a fit to [x°(T)]~ 3/2 . 



ture T c (x) is nonzero for all x < 1 and reaches zero only 
right at x = 1 (see Fig. [T]). Using optimal fluctuation 
theory, we have derived the functional forms of various 
thermodynamic observables in the tail of the smeared 
transition. We have also briefly discussed smeared classi- 
cal phase transitions that can occur in systems with cor- 
related disorder, and we have performed computer sim- 
ulations of a toy model that confirm and illustrate the 
theory. 

Although our results are qualitatively similar to those 
obtained for smeared phase transitions occurring at fixed 
randomness as a function of temperature or an appropri- 
ate quantum control parameter, the functional forms of 
observables are not identical. The most striking differ- 
ence can be found in the far tail of the transition. In the 
case of composition-tuning, the order parameter vanishes 
as a non-universal power of the distance from the end of 
the tail (x — 1), reflecting the fact that the minimum rare 
region size required for local magnetic order is finite. In 
contrast, if the transition occurs at fixed composition as 
a function of temperature or some quantum control pa- 
rameter, the order parameter vanishes exponentially, 2,4 
because the minimum size of an ordered rare region di- 
verges in the far tail. These differences illustrate the 
fact that the behavior of observables at a smeared phase 
transition is generally not universal in the sense of crit- 
ical phenomena; it depends on details of the disorder 
distribution and how the transition is tuned. Only the 
question of whether or not a particular phase transition is 
smeared is universal, i.e., determined only by symmetries 



and dimensionalities. 

Let us briefly comment on the relation of our theory to 
percolation ideas. The optimal fluctuation theory of Sec. 
IIIBI applies for compositions x larger than the percola- 
tion threshold of the A-atoms. Because the A-clusters 
are disconnected in this composition range, percolation 
of the A atoms does not play a role in forming the tail 
of the ordered phase at large x. Instead, distant rare 
regions are coupled via the fluctuations of the paramag- 
netic bulk phase and, in metallic magnets, via the RKKY 
interaction. Percolation does play a role, though, in the 
crossover between the inhomogeneous order in the tail of 
the transition and the bulk order at lower x. 

We note in passing that the behavior of a diluted sys- 
tem (where B represents a vacancy) with nearest-neighbor 
interactions is not described by our theory. In this case, 
the A-clusters are not coupled at all for compositions x 
larger than the A percolation threshold. Therefore they 
cannot align, and long-range order is impossible. As a 
result, the super-paramagnetic behavior of the locally 
ordered clusters extends all the way down to zero tem- 
perature. This was recently discussed in detail on the 
example of a diluted dissipative quantum Ising model. 21 

In the present paper, we have assumed that the A and 
B atoms are distributed independently over the lattice 
sites, i.e., we have assumed that there are no correlations 
between the atom positions. It is interesting to ask how 
the results change if this assumption is not fulfilled, for 
example because like atoms tend to cluster. As long as 
the correlations of the atom positions are short-ranged 
(corresponding to a finite, microscopic length scale for 
clustering), our results will not change qualitatively. All 
arguments in the optimal fluctuation theory still hold 
using a typical cluster of like atoms instead of a single 
atom as the basic unit. However, such clustering will 
lead to significant quantitative changes (i.e., changes in 
the non- universal constants in our results), as it greatly 
increases the probability of finding large locally ordered 
rare regions. We thus expect that clustering of like atoms 
will enhance the tail and move the phase boundary T c (x) 
towards larger x. A quantitative analysis of this effect 
requires explicit information about the type of correla- 
tions between the atom positions and is thus relegated 
to future work. 

Let us finally turn to experiment. Tails of the ordered 
phase have been observed at many quantum phase tran- 
sitions. However, it is often not clear whether these tails 
are an intrinsic effect or due to experimental difficul- 
ties such as macroscopic concentration gradients or other 
macroscopic sample inhomogeneities. Recent highly sen- 
sitive magneto-optical experiments on Sr1_2.Caa.RuO3 
have provided strong evidence for a smeared ferromag- 
netic quantum phase transition^ The behavior of the 
magnetization and critical temperature in the tail of the 
smeared transition agree well with the theory developed 
here. Moreover, the effects of clustering discussed above 
may explain the wide variation of the critical composition 
between about 0.5 and 1 reported in earlier studies^ - — 
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We expect that our smeared quantum phase transition no. DMR-0906566. 
scenario applies to a broad class of itinerant systems with 
quenched disorder. 
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